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ABSTRACT 

We employ measurements of the [a/Fe] ratio derived from low-resolution (R ~ 2000) spectra of 17,277 
G-type dwarfs from the SEGUE survey to separate them into likely thin- and thick-disk subsamples. Both 
subsamples exhibit strong gradients of orbital rotational velocity with metallicity, of opposite signs, -20 to 
-30 km s" 1 dex -1 for the thin-disk and +40 to +50 km s" 1 dex" 1 for the thick-disk population. The rotational 
velocity is uncorrelated with Galactocentric distance for the thin-disk subsample, and exhibits a small trend 
for the thick-disk subsample. The rotational velocity decreases with distance from the plane for both disk 
components, with similar slopes (-9.0 ± 1.0 km s" 1 kpc -1 ). Thick-disk stars exhibit a strong trend of orbital 
eccentricity with metallicity (about -0.2 dex" 1 ), while the eccentricity does not change with metallicity for 
the thin-disk subsample. The eccentricity is almost independent of Galactocentric radius for the thin-disk 
population, while a marginal gradient of the eccentricity with radius exists for the thick-disk population. Both 
subsamples possess similar positive gradients of eccentricity with distance from the Galactic plane. The shapes 
of the eccentricity distributions for the thin- and thick-disk populations are independent of distance from the 
plane, and include no significant numbers of stars with eccentricity above 0.6. Among several contemporary 
models of disk evolution we consider, radial migration appears to have played an important role in the evolution 
of the thin-disk population, but possibly less so for the thick disk, relative to the gas-rich merger or disk heating 
scenarios. We emphasize that more physically realistic models and simulations need to be constructed in order 
to carry out the detailed quantitative comparisons that our new data enable. 

Subject headings: Galaxy: disk — Galaxy: formation — Galaxy: kinematics and dynamics — Galaxy: structure 



1. INTRODUCTION 

The Milky Way's thick disk, first identified from fits of the 
vertical density profile of stars with a mix of exponential func- 
tions (Yoshii 1982; Gilmore & Reid 1983), differs in many 
ways from the thin disk, e.g., in its kinematics and chemical 
abundances. 

The scale height of the thick disk is about 1 kpc, while that 
of the thin disk is ~ 0.3 kpc. Typical thick-disk stars have gen- 
erally lower net orbital rotational velocities with larger veloc- 
ity dispersions (Majewski 1993; Chiba & Beers 2000; Robin 
et al. 2003; Soubiran et al. 2003; Parker et al. 2004; Wyse 
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et al. 2006), possess higher [a/Fe] ratioQ, and are older 
and more metal-poor than typical thin-disk stars (Bensby et 
al. 2003, 2005; Reddy et al. 2006, 2010; Fuhrmann 2008; 
Haywood 2008). 

Their higher [a/Fe] ratios and older ages imply that thick- 
disk stars were born earlier than most thin-disk stars, in an 
environment of rapid star formation, and that they have likely 
had more time to experience dynamical heating and secular 
processes such as scattering by perturbations in the disk. As a 
result of the multiple complex processes that thick-disk stars 
may have experienced during their lifetimes, consensus on the 
nature of the formation and evolution of the thick disk has yet 
to be reached. 

The currently discussed mechanisms for thick-disk forma- 
tion can be broadly divided into two groups - violent ori- 
gin and secular evolution. Among the models involving vi- 
olent origin, the heating scenario (e.g., Quinn et al. 1993; 
Kazantzidis et al. 2008) posits that the thick disk results from 
a pre-existing thin disk that has been dynamically heated by 
satellite mergers. In their simulations of this process, Villalo- 
bos & Helmi (2008) found that on the order of 10-20% of 
the stars in the thickened disk component were accreted from 
satellites, the rest being heated thin-disk stars. The accretion 
origin of the thick disk (e.g., Abadi et al. 2003) invokes the 
hypothesis that thick-disk stars were predominantly formed 
in dwarf-like galaxies, which were then directly assimilated 
into the thick disk from orbits that reached near the Galactic 
disk plane. Abadi et al. (2003) predicted that over 70% of 

12 The [a/Fe] ratio is often represented by an average of the [Mg/Fe], 
[Si/Fe], [Ca/Fe], and [Ti/Fe] ratios, which we adopt in this paper as well. 
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thick-disk stars were accreted from such disrupted galaxies. 
The third model among the violent origin class is that thick- 
disk stars may have formed in situ through chaotic mergers of 
gas-rich systems, prompting simultaneous early star forma- 
tion before and during the mergers (Brook et al. 2004, 2005, 
2007), and that thin-disk stars formed after the merger events 
settled down. 

Secular evolution by disk heating was first conceived by 
Spitzer & Schwarzschild (1953), who demonstrated that en- 
counters with molecular clouds could increase the velocity 
dispersion of late type, old stars. Barbanis & Woltjer (1967) 
also showed that spiral structures might be the cause of the 
larger velocity dispersion of older stars in the solar neigh- 
borhood. These ideas had been further developed by several 
studies (e.g., Fuchs 2001, and references therein). Although 
challenged by Jenkins (1992), disk heating by secular pro- 
cesses have recently regained attention, both observationally 
and theoretically, as possible thick-disk formation scenarios. 

Indeed, recent theoretical studies and simulations (Schon- 
rich & Binney 2009a, 2009b; Loebman et al. 2010) suggested 
that the thick disk might not require a violent origin, but rather 
could have formed by cumulative secular processes associated 
with the radial migration of stars. According to the migration 
theories (Sellwood & Binney 2002; Roskar et al. 2008a), stars 
in the Galactic disk can radially move from the inner (outer) 
to the outer (inner) regions due to resonant scattering by tran- 
sient spiral structure. Based on their simulations, Minchev 
& Famaey (2010) also suggested that long-lived spiral struc- 
tures, interacting with a central bar, could be responsible for 
the radial movements of stars in a disk galaxy. 

These proposed models predict various trends between 
the kinematic parameters and chemical abundances of disk- 
system stars, as well as between their kinematics and spatial 
distributions. For example, Schonrich & Binney (2008b) sug- 
gested that local, relatively metal-rich thin-disk stars, formed 
in the inner part of the disk and moved outward, while lo- 
cal metal-poor thin-disk stars were born in the outer disk and 
migrated inward to the solar radius, retaining information on 
the kinematic differences between the two populations. Thus, 
there should exist a gradient in the variation of rotational ve- 
locity with metallicity; evidence for such a behavior has been 
claimed observationally by Haywood (2008). Models of disk 
heating via satellite mergers (Villalobos et al. 2010) result in 
proposed relationships between rotational velocity and Galac- 
tocentric distance and distance from the Galactic plane. Gas- 
rich merger models (Brook et al. 2007) also predict a gradient 
of rotational velocity with Galactocentric radius for disk stars 
near the solar radius. 

Sales et al. (2009) proposed that the distribution of orbital 
eccentricities for nearby thick-disk stars could be used to pro- 
vide constraints on the various suggested formation models. 
A number of recent papers also have employed this frame- 
work to study possible origins of the thick disk, based on data 
from several large spectroscopic surveys. For example, Wil- 
son et al. (201 1) have explored data from the RAdial Velocity 
Experiment (RAVE; Steinmetz et al. 2006), while Dierickx 
et al. (2010) used data from the seventh public release of 
the Sloan Digital Sky Survey (SDSS DR7; York et al. 2000; 
Abazajian et al. 2009). The study of Casetti-Dinescu et al. 
(201 1) combined RAVE data with newly available proper mo- 
tions from the fourth release of the Southern Proper Motion 
Catalog (SPM4; Girard et al. 201 1). We discuss their analyses 
and conclusions further below. 

Most previous observational studies that have sought to test 



the various correlations predicted by the models mentioned 
above have used methods of assigning individual stars to 
membership in the thin- and thick-disk populations (based on 
a given star's location or kinematics) that introduce manifest 
biases that can confound interpretations (as previously noted 
by Schonrich & Binney 2009b and Loebman et al. 2010). As 
the chemical signatures of a star are substantially less vari- 
able properties than its spatial position or velocities over its 
lifetime, it is instead desirable to classify disk stars into their 
likely components according to their chemistry. 

Among the various chemical abundance ratios that might 
be explored for this purpose, the [a/Fe] ratios appear particu- 
larly useful. These ratios can be relatively easily measured 
(as described below), and have been proven to well sepa- 
rate thick-disk stars from thin-disk stars. It is known that, 
at least in the solar neighborhood (where essentially all pre- 
vious studies have been conducted), thick-disk stars are on 
average enhanced in their [a/Fe] ratios by +0.2 to +0.3 dex 
relative to their thin-disk counterparts at a given [Fe/H]. Lo- 
cal kinematically selected thin- and thick-disk samples based 
on probabilistic membership assignments have confirmed this 
enhancement of [a/Fe] (e.g., Bensby et al. 2003; Reddy et al. 
2006). Fuhrmann (1998, 2008) also demonstrated that dwarfs 
in his volume-limited sample could be clearly separated into 
two populations as a function of [Fe/H] - one associated with 
high [Mg/Fe] and the other with low [Mg/Fe] . The elemen- 
tal abundance patterns of the stars with low [Mg/Fe] ratios 
and high [Mg/Fe] ratios are very similar to the kinematically 
selected thin- and thick-disk samples. The recent study by 
Nissen & Schuster (2010) demonstrated that nearby dwarfs 
with halo kinematics could be separated into two groups on 
[a/Fe]. They proposed that the high-a stars may have been 
born in the disk or bulge of the milky way and heated to halo 
kinematics by merging satellite galaxies or else were simply 
members of the early generations of halo stars born during the 
collapse of a proto-Galactic gas cloud, while the low-a stars 
may have been accreted from dwarf galaxies. Clearly, [a/Fe] 
ratios also provide valuable information on the timescales and 
intensities of star formation in the populations involved. 

In this study we make use of the first set of [a/Fe] ratios ob- 
tained for a large sample of low-resolution (R ~ 2000) spectra 
from the Sloan Extension for Galactic Understanding and Ex- 
ploration (SEGUE; Yanny et al. 2009). As shown by Lee et 
al. (2011), for stars with SDSS/SEGUE spectra of signal-to- 
noise (S/N) ratios greater than 20 A" 1 , and with temperatures 
in the range 4500 K < T e ff < 7000 K, one can estimate [a/Fe] 
with an accuracy of better than 0. 1 dex (as derived by com- 
paring [a/Fe] estimates from the analysis of moderately high- 
resolution and medium-resolution spectra with those obtained 
from application of our methods to the same stars). This en- 
ables a chemical separation of the disk system into likely 
thin- and thick-disk populations. In this paper we explore the 
observed correlations of rotational velocity and orbital eccen- 
tricity with metallicity, Galactocentric distance and distance 
from the Galactic plane, as well as the orbital eccentricity dis- 
tributions for the individual populations, and compare with 
the predictions of the radial migration, gas-rich merger, ac- 
cretion, and dynamical heating models. Since we believe that 
direct quantitative comparisons with the predictions made by 
various models (or simulations) mentioned in this study are 
somewhat premature, we emphasize the more qualitative as- 
pects of these comparisons. 

This paper is outlined as follows. In Section 2 we present 
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the G-dwarf sample from SEGUE, describe various cuts im- 
posed on the sample to obtain a refined disk dwarf sample, 
and discuss the calculations used to derive their space mo- 
tions and orbital eccentricities. Section 3 describes how we 
assign membership of the stars into either the thin- or thick- 
disk populations. Results from our G-dwarf sample and dis- 
cussion of comparisons of our results with the predictions of 
various contemporary disk formation and evolution scenarios 
follow in Sections 4 and 5, respectively. A summary and our 
conclusions follow in Section 6. 

2. SELECTION OF LOCAL DWARF STARS 

2.1. The SEGUE G-dwarf Sample 

Our initial sample comprises low-resolution (R —2000) 
spectra of -63,000 stars from SDSS Data Release 8 (DR8; 
Aihara et al. 2011), obtained during the SEGUE sub-survey, 
which were originally targeted as G-dwarf candidates (with 
colors and magnitudes in the range 0.48 < (g—r)o < 0.55 and 
ro < 20.2). As a result of the simple sampling function, this 
dataset is expected to be relatively unbiased with respect to 
chemistry, and completely unbiased with respect to kinemat- 
ics. In order to obtain a subsample of disk stars with the most 
reliably estimated physical quantities we apply several addi- 
tional cuts. 

First, we exclude stars lacking information on their stellar 
parameters (effective temperature, T t g, surface gravity, log g, 
and metallicity, [Fe/H]), radial velocities, or proper motions. 
The stellar atmospheric parameters were determined by the 
most recent version of the SEGUE Stellar Parameter Pipeline 
(SSPP; Lee et al. 2008a, 2008b; Allende Prieto et al. 2008; 
Smolinski et al. 2011); typical external errors in these esti- 
mates are 180 K in T e s, 0.24 dex in log g, and 0.23 dex in 
[Fe/H] (Smolinski et al. 2011). It has been shown that shifts 
in the SSPP-derived estimates of [Fe/H] and [a/Fe] caused by 
the presence of unrecognized spectroscopic binaries are gen- 
erally small (Schlesinger et al. 2010). Although the typical 
uncertainty of the radial velocity varies with the S/N ratio of 
a spectrum, it is less than 5 km s" 1 for the great majority of 
stars in our sample. Proper motion information was obtained 
based on the procedures described by Munn et al. (2004); 
the systematic error noted by Munn et al. (2008) has been 
corrected (final typical errors are 3-4 mas yr" 1 ), and we have 
adopted the Munn et al. recommendations for maximum fit 
residuals and minimum numbers of epochs considered in or- 
der to obtain the most reliable proper motions. In this regard, 
see also Bond et al. (2010), who investigated the systematic 
errors in Munn et al. (2008) by comparison with the expected 
null proper motions of SDSS quasars. 

Distances to individual stars are estimated using a cali- 
brated set of stellar isochrones (An et al. 2009a), follow- 
ing the prescription in An et al. (2009b). After correcting 
photometry for dust extinction, main-sequence fitting is per- 
formed simultaneously on three different color-magnitude di- 
agrams (CMDs), with r as a luminosity index, and g — r, g-i, 
and g — z as color indices, respectively. We adopt an SSPP- 
derived [Fe/H] in the distance estimation for each star, and fix 
a stellar age and [a/Fe] of the model at a given [Fe/H], as- 
suming a linear relationship between [Fe/H] and these quan- 
tities (see An et al. 2009b). Distance estimates obtain using 
[a/Fe] from this assumption may not be internally consistent 
with analyses based on the SSPP-determined [a/Fe], but even 
a ~0.1 dex difference in [a/Fe] has a negligible impact on the 
derived distances (—0.01 mag in distance modulus). We also 



limit models in the fitting to log g > 4.2 to minimize possible 
distance bias from stellar age effects near the main-sequence 
turn-off. An inter-comparison of results from the three CMDs 
suggests that the internal error in the distance modulus is —0. 1 
mag; an additional ~0.1 mag error is expected from the com- 
bined errors in age, [Fe/H], [a/Fe], and E(B-V). This sug- 
gests that the associated distance-modulus error is ~0.14 mag 
for individual stars. The effects of binarity are more difficult 
to quantify, and are not included in this error estimate (see An 
et al. 2007; Sesar et al. 2008). 

The [a/Fe] ratio is derived following the procedures de- 
scribed by Lee et al. (2011). Briefly summarizing, Lee et 
al. first generated a grid of synthetic spectra, covering 4000 K 
< r eff < 8000 K in steps of 250 K, 0.0 < log g < 5.0 in steps 
of 0.2 dex, -4.0 < [Fe/H] < +0.4 in steps of 0.2 dex, and 
-0.1 < [a/Fe] < +0.6, in steps of 0.1 dex, then determined 
[a/Fe] by searching the grid for a synthetic spectrum that best 
matches a given SDSS/SEGUE spectrum (in regions that are 
most influenced by [a/Fe]). By comparing with a set of mod- 
erately high-resolution (R = 15,000) and medium-resolution 
(R = 6000) spectra of SDSS/SEGUE stars, they demonstrated 
the ability to measure [a/Fe] from SDSS/SEGUE spectra 
(with S/N > 20 A -1 ) with uncertainties less than 0.1 dex, for 
stars with atmospheric parameters in the range T e g = [4500, 
7000] K, log g = [1.5, 5.0], and [Fe/H] = [-1.4, +0.3], over the 
full range of [a/Fe] considered. For stars with [Fe/H] < -1 .4, 
slightly higher S/N was required to achieve this precision (S/N 
> 25 A" 1 ). 

In order to assemble a local dwarf sample, we only include 
stars with distances, d, less than 3 kpc from the Sun, and with 
log g > 4.2. These cuts ensure that we are selecting likely 
dwarfs from which we can obtain accurate space motions (i.e., 
that do not suffer from severe degradation due to propagation 
of proper motion errors at larger distances). In order to per- 
form a confident separation of the thin- and thick-disk popula- 
tions on the basis of [a/Fe], we further require that the spectra 
of the dwarf stars included in our analysis have S/N > 30 A" 1 . 
This conservative cut on S/N ensures not only high quality es- 
timates of [Fe/H] and [a/Fe], but also that our program stars 
have small errors in estimated radial velocity (less than 5 km 
s" 1 ). 

2.2. Calculations of Space Motions and Orbital Eccentricity 

With information on the distances, radial velocities, and 
proper motions for our program stars in hand, we then derive 
the U, V, W space velocity components. We apply (U,V,W)q 
= (11. 10,12.24,7.25) kms" 1 (Schonrichet al. 2010) to adjust 
for the solar peculiar motions with respect to the Local Stan- 
dard of Rest (LSR). For the purpose of our analysis, we also 
make use of the rotational velocity around the Galactic center 
in a cylindrical coordinate system, V«^, calculated assuming 
Rq = 8.0 kpc and Vlsr = 220 km s" 1 . The Galactocentric 
distance projected onto the Galactic plane, R, and the verti- 
cal distance from the Galactic plane, |Z|, are also obtained. 
In addition, by adoption of an analytic Stackel-type gravita- 
tional potential (which includes a flattened, oblate disk and 
a spherically-shaped massive dark halo; see Chiba & Beers 
2000), we compute r apo (r per i), the maximum (minimum) dis- 
tance from the Galactic center that a star reaches during its or- 
bit, as well as the orbital eccentricity, e, defined as (r apo -r per i) 
AVapo + Tperi)- Errors in the derived kinematics and orbital pa- 
rameters for each star due to propagation of the errors in the 
observed quantities (mostly dominated by distance and proper 
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Figure 1. Distributions of metallic! ties, [Fe/H] (top left), rotational velocities, V^, (top right), distances from the Galactic plane, |Z (bottom left), and Galacto- 
centric distances projected onto the plane, R (bottom right) for the final selected sample. The solid line indicates our final full dwarf sample, while the dashed red 
and dotted-dash blue lines are the thin- and thick-disk subsamples, respectively, with stars assigned by the procedures described in Section 3. 



motion errors) are determined by 1000 realizations of a Monte 
Carlo simulation. 

We next remove stars from our sample with derived rota- 
tional velocities relative to the Galactic center less than V<j> = 
+50kms _1 , with [Fe/H] <— 1.2, and located outside the range 
7 < R < 10 kpc, in order to minimize contamination from the 
halo and outer-disk components. 

Finally, we perform a simple check on the likely remaining 
halo contamination in our sample following the prescription 
of Bensby et al. (2003). For calculation of the approximate 
disk and halo star fractions (assuming our data is represen- 
tative of the local solar neighborhood densities), we adopt 
the local stellar densities, velocity dispersions in U, V, and 
W, and the asymmetric drifts listed in their Table 1, assum- 
ing the space velocities of the thin-disk, thick-disk, and halo 
stars are distributed as Gaussians. Based on these probabil- 
ity distributions, we reject stars that have greater likelihood of 
belonging to the halo than to the disk system. This check re- 
moves only about 59 additional stars from the sample, show- 
ing that the above selection criteria for thin- and thick-disk 
stars are quite reasonable. We also experimented with the ap- 
plication of slightly different scale heights for describing the 
variation of halo stellar densities with |Z|, but the above re- 
sult appears quite robust. Note, however, that these various 
cuts do not necessarily eliminate contamination by members 
of the so-called metal-weak thick-disk (MWTD) population, 
which Carollo et al. (2010) have shown exhibits metallicities 
in the range -1.7 < [Fe/H] < -0.7, and a prograde rotation of 
Vcf, ~ +100 to +150 km s" 1 . We comment on any evidence for 
MWTD contamination below. 



Summarizing the criteria used for our sample selection, sur- 
viving program stars satisfy d < 3 kpc, log g > 4.2, S/N > 30 
A" 1 , V > +50 km s _1 , [Fe/H] > -1.2, 7 < R < 10 kpc, and 
possess greater probability of belonging to the disk system 
than to the halo. The surviving sample from the above cuts 
numbers ^17,300 stars. Figure Q] shows the distributions of 
[Fe/H], Vrf,, \Z\, and R for the final dwarf sample (solid lines), 
before and after further division based on the derived [a/Fe] 
ratios into the thin- and thick-disk populations, as described 
below. Note that in the bottom left panel of Figure Q]has only 
a small number of stars with |Z| < 0.2 kpc, owing to the bright 
limit of SDSS imaging (g > 14.0). Thus, our analysis in the 
following sections may be only valid for the thin-disk popu- 
lation with |Z| > 0.2 kpc, rather than the young(er) thin disk 
closer to the plane. 

3. DIVISION OF THE SAMPLE ON [a/Fe] INTO THIN- AND 
THICK-DISK POPULATIONS 

3.1. Dividing Scheme and its Efficiency 

As mentioned previously, since a stellar population's kine- 
matics and spatial distributions can be modified over time 
(especially in the disk system), while a (dwarf) star's atmo- 
spheric chemical abundance is essentially invariant (except in 
unusual circumstances, such as binary mass transfer from an 
evolved companion), we make use of the estimated [a/Fe] ra- 
tio as a reference to separate the thin- and thick-disk popula- 
tions. This choice is also strongly motivated by the apparent 
bi-modal distribution of stars in the [a/Fe] and [Fe/H] plane 
seen in Figure [2] 

For the purpose of the present analysis, our dwarf sample 
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Figure 2. Distribution of logarithmic number densities, in the [a/Fe] vs. 
[Fe/H] plane, over-plotted with equidensity contours. Each bin is 0.025 dex in 
[a/Fe] by 0.05 dex in [Fe/H], and is occupied by a minimum of 20 stars. The 
median occupancy is 70 stars. The solid line is the fiducial for division into 
likely thin- and thick-disk populations; the dashed lines located ±0.05 dex in 
[a/Fe] on either side of the solid line indicate the adopted dividing points for 
the high-[a/Fe] (upper-dashed) and low-[a/Fe] (lower-dashed) stars in our 
sample. 

is split into likely thin-disk (with low [a/Fe]) and thick-disk 
(with high [a/Fe]) populations, based on the following 
scheme: 

I) For stars with [Fe/H] > -0.8 

• thin disk, if [a/Fe] < -0.08[Fe/H] + 0.15 

• thick disk, if [a/Fe] > -0.08-[Fe/H] + 0.25 

II) For stars with [Fe/H] < -0.8 

• thin disk, if [a/Fe] < +0.214 

• thick disk, if [a/Fe] > +0.314 

This division into the thin- and thick-disk populations is de- 
vised based on examination of the distribution of number den- 
sities in the [a/Fe] vs. [Fe/H] plane, shown in Figure [2] Note 
how well the populations appear to separate above and below 
the solid line in this figure, which is our adopted fiducial. 

The dashed lines located ±0.05 dex in [a/Fe] above and 
below the fiducial solid line in Figure [2] indicate the dividing 
points for the high-[a/Fe] (thick-disk) and low-[a/Fe] (thin- 
disk) stars. Note that this leaves a gap of 0.1 dex in [a/Fe] 
between the thin- and thick-disk dividing lines. This choice 
serves to reduce the number of misclassified stars that may 
arise from observational errors in their measured [a/Fe]. The 
dashed red line in Figure [Tj shows the thin-disk subsample, 
whereas the dotted-dash blue line is for the thick-disk sub- 
sample, classified by the dividing schemes described above. 
From this figure, one can roughly read off the ranges and peak 
values of the estimated and derived parameters for each sub- 
sample. 

To check on the efficacy of the chemical separation of the 
disk populations through use of the [a/Fe] ratio, we have in- 
vestigated the variation of the U, V, W velocity dispersions 
of our sample with [a/Fe]. It is well known that the disper- 
sion of each velocity component increases with distance from 
the Galactic plane (as well as on age, on average). In any 
event, the thick-disk population exhibits substantially higher 
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Figure 3. [/, V, W velocity dispersions, as a function of [a/Fe], are shown 
from top to bottom, respectively. Well-defined trends of the velocity disper- 
sions can be seen, providing kinematic confirmation that the proposed separa- 
tion of the thin-disk and thick-disk populations on the basis of [a/Fe] works 
very well. Each dot contains 200 stars, and the error bar on each point is 
calculated by resampling these 200 stars with replacement 1 000 times. The 
vertical dashed lines provide references at [a/Fe] = +0.2 and +0.3, which 
roughly correspond to the divisions shown in Figure [2] 

dispersions than the thin-disk counterpart. Figure [3]shows the 
derived velocity dispersions of our sample as a function of 
[a/Fe]. It is readily apparent that, up to around [a/Fe] = +0.2, 
the dispersion of each velocity component increases moder- 
ately. Above [a/Fe] = +0.2 the gradients of the velocity dis- 
persions with [a/Fe] become somewhat steeper. Above [a/Fe] 
= +0.3, the magnitude of each velocity dispersion is larger by 
about 10 km s" 1 than for [a/Fe] < +0.2. As our thin-disk stars 
mostly have [a/Fe] < +0.2 and thick-disk stars possess [a/Fe] 
> +0.3, Figure [3] kinematic ally confirms that the division by 
[a/Fe] into the thin- and thick-disk populations is quite robust. 

3.2. Impact of Potential Metallicity Bias 

One may be concerned about biases in our initial sample 
which might arise from our (g - r)o target selection due to 
small, but non-zero, metallicity sensitivity of the stellar col- 
ors. Haywood (2001), for example, has previously suggested 
that targeting stars on the basis of their spectral types (e.g., G- 
dwarfs) preferentially selected more metal-poor than metal- 
rich stars, resulting in a metallicity distribution shifted lower 
by about 0. 1-0.3 dex. 

If our selected sample strongly favors metal-poor over 
metal-rich stars, this bias might produce misleading correla- 
tions between the parameters we are seeking to understand. 
For example, at least for the thick-disk population, previous 
studies have indicated that the observed stellar orbital rota- 
tional velocity decreases with declining metallicity. Thus, if 
biases have increased the relative numbers of metal-poor stars 
in the thick-disk subsample, the overall distribution of will 
be shifted to lower rotational velocity. However, it should be 
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Figure 4. Distributions of metallicities in different regions of distance from the Galactic plane (left panel) and the Galactic center (right panel) for the thin-disk 
(solid line) and thick-disk (dashed line) subsamples. For the thin-disk subsample, as the |Z| distance increases, the relative numbers of metal-poor stars slightly 
increases (in particular for the region |Z| > 1.5 kpc). For the thick-disk subsample, the metallicity distribution does not change with height above or below the 
Galactic midplane. The metallicity distribution mostly stays the same throughout the different cuts of Galactocentric distance for both the thin- and thick-disk 
subsamples. The bottom panels show the full samples of thin- and thick-disk stars, over the range 0.1 < |Z| < 3.0 kpc. Note that the y-axis scale is very different 
for each panel. The vertical dash-dotted lines are reference lines at [Fe/H] = 0.0 and [Fe/H] = -0.6. 



kept in mind that, because our sample does not suffer from 
kinematic bias, any correlations that we are seeking between 
kinematics and chemical abundances will not be affected by 
any potential metallicity bias, as long as the correlations are 
derived from ranges of R and |Z| that are sufficiently small 
that the correlations remain roughly constant over the regions 
considered. Furthermore, any kinematic trends with R and |Z| 
will not be affected by metallicity bias, as long as the metal- 
licity distributions of stars in different ranges of the spatial 
cuts do not vary significantly. 

Hence, instead of correcting for possible selection bias in 
our sample (which itself is a complex and tricky business), we 
instead seek to demonstrate that any presumed metallicity bias 
does not greatly impact the kinematic correlations with spatial 
parameters that we derive by examination of the metallicity 
distribution functions (MDFs) for both the thin- and thick- 
disk subsamples in various regions of R and |Z|. If there is 
any sudden change in the MDFs of each subsample between 
neighboring regions, that will be a sign that the subsample 
may not be suitable for deriving meaningful correlations be- 
tween kinematics and distance. 

For the thick disk, there is some existing evidence for the 
lack of a metallicity gradient with distance above the Galactic 
plane (e.g., Gilmore et al. 1995), or at most for only a small 
one, on the order of 0. 1-0.2 dex kpc" 1 (Ivezic et al. 2008). So, 
the shape of the observed MDFs at different heights should 
remain roughly constant. The left panel of Figure |4]displays 
the observed MDFs for both the thin- and thick-disk subsam- 
ples in different bins of |Z| distance. From inspection, the 



relative numbers of metal-rich stars in the thick-disk subsam- 
ple do not grossly change with different cuts on height above 
the plane, so false kinematic trends with |Z| are not expected 
to arise. Quantitatively, the fraction of the stars with [Fe/H] 
< -0.6 for the thick-disk subsample is 0.51, 0.51, 0.58, and 
0.64 from the first to the fourth panel (a resulting metallic- 
ity gradient with |Z| of -0.029 ± 0.019 dex kpc" 1 ), consistent 
with the expectation from previous work. These results are 
also in line with the findings of Haywood (2001), who noted 
that the metallicity bias arising from the color selection was 
significant for stars with [Fe/H] > -0.2, but not for relatively 
metal-poor thick-disk stars. 

However, this seems not to be the case for the thin-disk sub- 
sample shown in the left panel of Figure [4] At heights above 
|Z| = 1.5 kpc (already many thin-disk scale heights above the 
plane), we notice that more metal-rich stars have dropped out 
of the distribution, compared to the upper three panels. This 
is quantitatively confirmed by examination of the fraction of 
stars with [Fe/H] < -0.2 (0.45, 0.46, 0.58, and 0.71 from the 
first to the fourth panel). This may be a natural consequence 
of selecting the sample without consideration of the |Z| dis- 
tance, since at greater heights thick-disk stars are expected 
to dominate. In other words, some of the stars in the metal- 
poor tail of the thin-disk subsample may in reality belong to 
the thick disk, but they have been misclassified as thin-disk 
stars due to errors in the estimated [a/Fe]. Indeed, consider- 
ing the distribution of [a/Fe] for the stars with [Fe/H] < -0.3 
and |Z| > 1.5 in the thin-disk subsample, many of the stars 
have [a/Fe] > +0.15. Thus, moderately a-enhanced stars at 
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this height may mostly belong to the thick disk rather than to 
the thin disk, with a much lower probability of old thin-disk 
membership. 

Simple experiments support the above argument. Accord- 
ing to Lee et al. (201 1), the error in [a/Fe] at S/N = 30 A" 1 is 
about 0.08 dex. If we assume this is a reasonable estimate of 
a la random error in [a/Fe] for all stars with |Z| > 1.5 kpc, 
and simultaneously perturb the measured [a/Fe] and [Fe/H] 
values according to a normal distribution with a = 0.08 dex 
and the measured uncertainty in [Fe/H], the total number of 
stars classified as members of the thin-disk component falls 
to about 85, after averaging the results of 1000 different real- 
izations. This is substantially smaller than the 201 stars that 
are claimed to be present. Thus, it is valid (within statistical 
fluctuations) to say that the low- [a/Fe] metal -poor thin-disk 
stars in this |Z| distance region are likely spurious, and are 
found at roughly the expected level of contamination. More- 
over, as the total number of thin-disk stars in this most distant 
region is rather small, we expect the impact of such stars on 
our analysis to be minimal. 

Concerning the variation of the MDFs for cuts in the R dis- 
tance, many previous studies (e.g., Nordstrom et al. 2004; 
Holmberg et al. 2007; Andrievsky et al. 2004; Lemasle et 
al. 2007; Sestito et al. 2008; Magrini et al. 2009) that used 
open clusters and/or field stars to derive radial metallicity gra- 
dients report rather small gradients, on the order of -0.05 to 
-0.1 dex kpc" 1 . Thus, we expect the MDFs in our sample 
to appear similar in different bins of radial distance from the 
Galactic center. The right panel of Figure |4]shows the MDFs 
of the thin- and thick-disk subsamples for several slices in R. 
From inspection, the overall shapes of the MDFs remain the 
same through the fourth panel, although there is a small drop 
in the numbers of thin-disk stars around [Fe/H] ~ -0.1 in the 
distance range 9.0 < R < 10.0 kpc. Calculating the fraction 
of the stars with [Fe/H] < -0.6 for the thick-disk subsample, 
we obtain 0.56, 0.54, 0.58, and 0.61 from the first to the fourth 
panel, a resulting radial metallicity gradient of -0.008 ±0.004 
dex kpc" 1 . For the thin-disk subsample, the fraction of the 
stars with [Fe/H] < -0.2 is 0.42, 0.47, 0.48, and 0.57 from 
the first to the fourth panel. The derived radial metallicity 
gradient is -0.043 ± 0.004 dex kpc" 1 ; these slopes do not dif- 
fer significantly from the previous studies. Thus, the results 
of these two exercises suggest that it is highly unlikely that 
any potential metallicity bias in our sample can greatly affect 
our derived kinematic correlations with R and |Z|. 

Potential biases also might depend on the age distribution 
of our sample, which is not known at present, and in any case 
is difficult to quantify. As the narrow color range applied to 
originally select the G-dwarfs for spectroscopic follow-up in 
SEGUE also preferentially selects certain age ranges on the 
main sequence, we might expect that this bias might con- 
tribute at some level to the observed trends (e.g., rotational 
velocity versus metallicity) that we are seeking to understand. 
However, as noted by Haywood (2001), age bias is expected 
to be even less important than the metallicity bias that we have 
already shown to have minimal effect. 



4. RESULTS OF THE OBSERVATIONS 

In this section we use our local G-dwarf sample to examine 
the observed gradients of with [Fe/H], R, and |Z|, as well 
as trends of e with [Fe/H], R, and |Z| for the thin-disk and 
thick-disk populations as identified above. 



140 190 240 




Figure 5. Distribution of rotational velocities (V^>, top panel) and the average 
orbital radii (Rmem, bottom panel) for our G-dwarf sample in the [a/Fe] vs. 
[Fe/H] plane. As in Figure [2] the dividing lines for the thin- and thick-disk 
subsamples are shown. Each bin has a size of 0.025 dex in [a/Fe] by 0.05 dex 
in [Fe/H], and is occupied by a minimum of 20 stars. The median occupancy 
is 70 stars. Each bin represents a 3<r-clipped mean of so that outliers in 
each bin do not significantly affect the average behavior. 



4. 1 . Correlations between Rotational Velocity and Metallicity 

The top panel of Figure [5] shows a color-coded distribution 
of in the [a/Fe] vs. [Fe/H] plane for our G-dwarf sam- 
ple. Detailed examination of this panel (as well as Figure |2]i 
reveals a metal-poor tail for the low-[a/Fe] stars (< +0.2), 
which we associate with the thin disk, extending down to 
[Fe/H] = -0.7. This already implies that the thin disk may 
not be well-described by a single metal-rich population with 
a peak around [Fe/H] = -0.2. We also notice from this panel 
that a higher rotational velocity is observed in the region of 
the metal-poor thin disk ([a/Fe] < +0.2 and [Fe/H] < -0.3), 
suggesting a negative trend of with [Fe/H]. In contrast, 
the high-[a/Fe] stars (> +0.3, which we associate with the 
thick disk) apparently exhibit a strong positive trend of 
with [Fe/H]. We investigate these trends quantitatively below. 

The bottom panel of Figure [5] displays the distribution of 
mean orbital radii (/? mea n) in the [a/Fe] vs. [Fe/H] plane. It is 
clear that the stars we associate with the thick-disk population 
exhibit smaller mean orbital radii than those with associated 
with the thin disk. In addition, the metal-poor thin-disk stars 
possess larger mean orbital radii than the dominant metal-rich 
thin-disk stars. 

Looking at the results from other observational work, a re- 
cent study by Navarro et al. (2010) obtained a slightly differ- 
ent result for their thin-disk subsample. These authors found 
little or no correlation between V$ and [Fe/H] for their thin- 
disk stars (defined by [Fe/H] > -0.7 and [a/Fe] < +0.2), al- 
though the subset of their thin-disk subsample with available 
Eu abundances (so that potential thick-disk or halo stars could 
be rejected) exhibits a very similar pattern to that we identify 
here. One should also keep in mind the possibility of effects 
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Figure 6. Distribution of rotational velocities for the low-a, thin-disk sub- 
sample (top panel) and the high-a, thick-disk subsample (bottom panel) in 
the |Z| vs. [Fe/H] plane. As in Figure[3] each bin has a size of 0.025 dex in 
[a/Fe] by 0.05 dex in [Fe/H], and is occupied by a minimum of 20 stars. The 
median occupancy is 47 stars. Each bin represents a 3<r-clipped mean of Vj, 
so that outliers in each bin do not significantly affect the average behavior. 



from selection biases in their sample, as it was based on an 
assembly of stars that included kinematically selected targets. 

Haywood (2008) separated thin-disk stars with [Mg/Fe] 
< +0.2 from thick-disk stars with [Mg/Fe] > +0.2 in the spec- 
troscopic sample of Soubiran & Girard (2005), and found 
an increasing trend of the mean orbital radii with decreas- 
ing metallicity for the thin-disk population, along with a de- 
creasing trend of mean radii with decreasing metallicity for 
the thick-disk population (see his Figure 3). He claimed that 
this tendency resulted from stars that migrated from the inner 
and outer disk. These behaviors are qualitatively in very good 
agreement with our findings in the bottom panel of Figure [5] 

Rocha-Pinto et al. (2006) also reported a similar behavior 
between the mean orbital radii and the chemical abundances 
in their volume-complete sample of 325 late-type dwarfs. 
Their primary results were that, as the difference in the dis- 
tance between the mean orbital radius and the solar radius 
increases, the abundances of Fe, Na, Si, Ca, Ni, and Ba all 
decrease. This relationship between the chemical abundances 
and the mean orbital radii could be accounted for by radial 
displacements of the stars involved. 

It is quite remarkable that all of the observed behaviors of 
the mean orbital radii from our G-dwarf sample agree so well 
with several previous observational studies (based on much 
smaller samples). 

Figure [6] indicates that there exists a clear gradient of 
with [Fe/H] at any given |Z| distance, for both the thin-disk 
subsample (top panel) and the thick-disk subsample (bottom 
panel). Figure [7] displays the observed gradients at different 
heights above the plane for both subsamples. Similar slopes 
of V<j, for both the low- and high- [a/Fe] stars are obtained 
for the various slices in |Z| distance, although the slope of 
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•Nm, = 1947, AV,/A[Fe/H] = -22.2 ± 2.7 0.1 < IZI < 0.5 kpc 
□N UcII = 465, AV,/A[Fe/H] = +43. 1 ±8.8 




•Nm, =4829, AV,/A[Fe/H] = -23.5 ±1.9 0.5 < IZI < 1.0 kpc 
□N skk = 2324, AV/AfFe/H] = +44.0 ± 4. 1 




• N lli0 = 1085, AV„/A[Fe/H]= -31.2 ±5.2 1.0 < IZI < 1.5 kpc 
□N Hcl = 1803, AV/AfFe/H] = +47.1 ± 5.1 




•N 1110 = 201, AV„/A[Fe/H] = -48.4 ± 12.9 1.5 < IZI < 3.0 kpc 
n N UA = 994, AV,/A[Fe/H] = +37.1 + 8.1 
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Figure 7. Rotational velocity gradients with metallicity for different slices 
in distance from the Galactic plane, for stars assigned to the thin-disk (black 
dots) and thick-disk (open squares) populations. Each dot represents a 3<r- 
clipped average of 100 stars. Note that the Icr error bars associated with 
each point are very small (on the order of 3-5 km s -1 ), so for visualization 
purposes 2cr error bars are plotted instead. The error bars are calculated by 
resampling the 100 stars in each (with replacement) 1000 times. The bottom 
panel shows the results for the full samples of stars considered. Note that 
although binned data are shown for clarity, estimates of the slopes and their 
errors are obtained for the full unbinned data (see text). 

the thick-disk subsample becomes shallower at larger distance 
(fourth panel), and slightly steeper for the thin-disk subsam- 
ple (which only includes 201 stars). The slopes are obtained 
by performing least-squares fits to the unbinned samples of 
each population; the uncertainties in the slopes are calculated 
by resampling each population with replacement 1000 times, 
in that process, the [Fe/H] and V$ values are perturbed simul- 
taneously by draws from a normal distribution using the 1 a 
errors in each quantity. The data shown in the kinematic cor- 
relation plots are binned only for clarity. 

If significant contamination of our thick-disk subsample by 
unrecognized MWTD stars were present, we might expect the 
slope of the correlation of with [Fe/H] to increase with 
distance above the plane, due to the greater velocity lag, larger 
scale height, and lower metallities of the MWTD component 
compared with that of the canonical thick disk (Carollo et al. 
2010). That is, at larger distances from the plane and at lower 
[Fe/H], the mean would be expected to be lower than it 
would be for a pristine thick-disk sample. We see no evidence 
for steepening of the gradient in Figure [7] Note that this is not 
to be taken as a contradiction with the Carollo et al. (2010) 
results, as those considered a different sample of stars, most of 
which were substantially lower metallicity than considered in 
the present study, and explored larger heights above the plane. 
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A gradient of about -20 to -30 km s 1 dex on average, is 
shown to exist for the thin-disk subsample, and a strong gradi- 
ent of +40 to +50 km s" 1 dex -1 for the thick-disk subsample. 
This result for the thick-disk population agrees with the claim 
of Spagna et al. (2010), who derived a similar slope using F-, 
G-, and K-type dwarfs from SDSS DR7. However, this find- 
ing clearly contradicts the results of Ivezic et al (2008), who 
found little correlation between and [Fe/H]. We discuss a 
possible resolution to this discrepancy in the Appendix. 

In order to check how uncertainties in the parameters, V^, 
[Fe/H], [a/Fe], and |Z| affect our derived gradients of over 
|Z| for both the thin- and thick-disk subsamples, we have per- 
formed a simple Monte Carlo experiment. Assuming a nor- 
mal distribution with width set by the la estimated error for 
each parameter, random changes in each parameter were ap- 
plied to 1000 realizations of each subsample (over the full 
range in |Z|). We obtained average gradients of -19.5 ±1.0 
km s" 1 dex" 1 for the thin-disk subsample and +43.4 ± 1 .8 km 
s" 1 dex" 1 for the thick-disk subsample, in good agreement 
with those shown in the bottom panel of Figure [7] We also 
have not found any notable correlations among the uncertain- 
ties in the parameters, which could possibly affect the derived 
gradients. Thus, we believe that our derived gradients of 
over [Fe/H] are not grossly affected by errors in the derived 
parameters. 

4.2. Rotational Velocity Gradients with Distance from the 
Galactic Center and Galactic Plane 

Figure[8]shows the overall trends of rotational velocity with 
distance from the Galactic center (top panel) and with vertical 
distance from the plane (bottom panel) for the thin-disk (black 
dots) and thick-disk (open squares) populations. Inspection 
of the top panel of this figure indicates only a negligible ro- 
tational velocity gradient for the thin-disk subsample (only 
—0.1 ± 0.6 km s" 1 kpc" 1 ), consistent with a flat rotation curve 
in the solar neighborhood. The asymmetric drift is about 10 
km s" 1 at the solar radius, as found by previous work (e.g., 
Soubiran et al. 2003). A small gradient of -5.6 ± 1.1 km 
s kpc" 1 is found for the thick-disk subsample, which lags 
the Vlsr by —40 km s" 1 , not far from the lag of 5 1 km s" 1 ob- 
tained by Soubiran et al. (2003). Note that even if we include 
in the analysis the stars with < V$ < 50 km s" 1 that were 
eliminated in our original selection, we obtain very similar 
asymmetric drifts and gradients as for the case of a sample of 
stars that does not include them. The slopes and their uncer- 
tainties are obtained by the same way as for the gradient of 
with [Fe/H]. 

The bottom panel of Figure|8]shows that the gradients of 
with |Z| distance are very similar (about -9.0 km s" 1 kpc" 1 ) 
for both the thin- and thick-disk subsamples. The difference 
in V$ (the velocity lag) for the high- [a/Fe] stars relative to 
the low- [a/Fe] stars is almost constant, —30 km s" 1 at any 
given |Z| distance. This again suggests that contamination 
from MWTD stars is not a major issue for our thick-disk sub- 
sample. 

Comparing with other recent studies, the vertical gradient 
of Vff, with |Z| for our thick-disk subsample, -9.4 ± 1.3 km 
s" 1 kpc" 1 , is smaller than that obtained by Casetti-Dinescu et 
al. (2011), -25 km s" 1 kpc" 1 , based on -4400 red clump 
metal-rich thick-disk stars covering the metallicity range -0.6 
< [Fe/H] < +0.5, that of Ivezic et al. (2008) from their SDSS 
sample (-29 km s" 1 kpc" 1 ), that of Girard et al. (2006), who 
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Figure 8. Rotational velocity gradients with Galactocentric radius (top 
panel) and with height above the Galactic plane (bottom panel), for the thin- 
disk (black dots) and thick-disk (open squares) subsamples. As in Figure[7j 
each dot represents a 3cr-clipped average of 100 stars, with error bars cal- 
culated by resampling 100 stars (with replacement) 1000 times. Note that 
although binned data are shown for clarity, estimates of the slopes and their 
errors are obtained for the full unbinned data (see text). 

derived a gradient of -30 km s" 1 kpc" 1 from a sample of about 
1200 red giants located in the range |Z| = 1-4 kpc, and as ob- 
tained by Chiba & Beers (2000; -30 km s" 1 kpc" 1 ) for the 
subset of their non-kinematically selected stars in the metal- 
licity range -0.8 < [Fe/H] < -0.6 within 2 kpc of the Galactic 
plane. Even if we cut our thick-disk subsample to include only 
stars with [Fe/H] > -0.6, we obtain a slope of -7.6 ± 1 .7 km 
s" 1 kpc" 1 , consistent, within 3<r, with that derived from the 
subsample without a metallicity restriction. 

It is interesting to note that, if we consider our entire thin- 
and thick-disk subsamples with |Z| > 1.0 kpc together, we 
find a vertical gradient of -19.4 ± 1.8 km s" 1 kpc" 1 , in bet- 
ter agreement with the previous studies. That is, the derived 
vertical gradient of the rotational velocity becomes substan- 
tially steeper when the stars are not divided according to their 
[a/Fe] ratios. We conclude that accurate determination of the 
vertical gradient of with |Z| for the thick disk requires ap- 
plication of a chemical separation criterion (other than simply 
[Fe/H]) to isolate the various components. 

Application of our simple Monte Carlo experiment with 
1000 realizations of the subsamples yielded average radial 
gradients of with R of +0.1 ±0.4 km s" 1 kpc" 1 and 
—3.7 ±0.7 km s" 1 kpc" 1 , and vertical gradients with \Z\ of 
-9.2 ±0.9 km s" 1 kpc" 1 and -8.2 ±0.9 km s" 1 kpc" 1 for 
the thin- and thick-disk subsamples, respectively, without any 
notable covariances between the errors in the parameters in- 
volved. These values are very close to those listed in Figure 

m 

4.3. Correlations of Stellar Orbital Eccentricities with 
Metallicity, Distance from the Galactic Center, and 
Height Above the Galactic Plane 

Figure|9]shows trends of orbital eccentricities (e) for the G- 
dwarf sample, as a function of [Fe/H], R, and |Z|, from the top 
to bottom panel, respectively. The black dots denote our thin- 
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Figure 9. Trends of eccentricities as a function of [Fe/H], R, and |Z|, from 
the top to bottom panels, respectively, for the thin-disk (black dots) and 
thick-disk (open squares) subsamples. As in Figure [7] each dot represents 
a 3<r-clipped average of 100 stars, with error bars calculated by resampling 
100 stars (with replacement) 1000 times. Note that although binned data are 
shown for clarity, estimates of the slopes and their errors are obtained for the 
full unbinned data (see text). The calculated gradients of e and their uncer- 
tainties are listed in Table[T] 

disk subsample, while the open squares indicate the thick-disk 
subsample. 

One outstanding feature from inspection of the three panels 
is that the overall distribution of the orbital eccentricities for 
the thick-disk stars is easily separable from that for the thin- 
disk population. The top panel suggests that the trend of the 
eccentricities for the thin-disk stars is independent of metal- 
licity, i.e., an almost fiat trend of e with [Fe/H], indicative of 
a very narrow distribution of eccentricity, with a peak around 
e ~ 0.14. On the other hand, the trend of e for the thick-disk 
subsample generally increases as the metallicity decreases. A 
slope of-0.192±0.010dex _I is obtained from a least-squares 
fit to the averaged points. 

The second panel also shows several interesting features. 
As in the top panel, there is not much correlation between 
e and R for the thin-disk subsample, although the behavior 
trends slightly higher below R = 7.5 kpc. The thick-disk stars 
generally exhibit an increasing trend of e with increasing R. 
The eccentricity distributions for the thin- and thick-disk pop- 
ulations merge at R ~ 7.0 kpc. 

The bottom panel shows that the eccentricities for both low- 
and high-[a/Fe] stars increase on average the farther away 
they are from the Galactic plane. In addition, similar to Figure 
[8] it is also noticed that the difference (about 0. 1 ) in the eccen- 
tricity between the low-[a/Fe] and high-[a/Fe] subsamples is 



constant at any given |Z|. 

A simple Monte Carlo experiment with 1000 realizations of 
the subsamples also reveals that the derived trends of the ec- 
centricities with [Fe/H], R, and |Z| above are not strongly af- 
fected by errors in the parameters involved, as the computed 
gradients of the eccentricities are within 3<r from those val- 
ues listed in Table [T] which quantitatively summarizes various 
correlations discussed in this section for the two subsamples. 

5. QUALITATIVE COMPARISONS WITH PREDICTIONS OF 
CONTEMPORARY MODELS OF DISK FORMATION 

It may be unwise to rely too strongly on the present pre- 
dictions of the suggested thick-disk formation models. This 
follows because, even though they are able to reproduce some 
aspects of the Milky Way's disk system, the predicted proper- 
ties are limited by large uncertainties with their treatment of 
star formation, the dynamical interaction of presumed satel- 
lites with the disk, unavoidable numerical effects, and the 
myriad set of assumptions that are required in their construc- 
tion. Thus, in this section, we compare our observational find- 
ings only with qualitative expectations from the published ra- 
dial migration, gas-rich merger, accretion, and disk heating 
models. It is our expectation that, as the models and simu- 
lations improve, these comparisons will increasingly be able 
to discriminate between the relative importance of the various 
formation scenarios. 

5.1. Correlations between Rotational Velocity and Metallicity 

According to the radial migration models (Sellwood & Bin- 
ney 2002; Roskar et al. 2008a; Schonrich & Binney 2009a; 
Minchev & Famaey 2010), the (presumably) metal -poor stars 
of the thin disk (which includes young, low-[a/Fe] stars) that 
were born in the outer disk move inward to the solar neigh- 
borhood, while the (presumably) metal-rich stars that formed 
in the inner disk migrate outward into the solar neighborhood 
(as the inner region of the disk has a higher stellar and gas 
density, and is rapidly chemically enriched, most of the stars 
should be metal rich). 

Schonrich & Binney (2009a) suggested that this radial 
movement can occur by two mechanisms: "blurring" and 
"churning". Blurring refers to the increase of eccentricities 
over time at a similar angular momentum due to scattering, 
e.g., on giant molecular clouds. Churning is mostly trig- 
gered by resonant scattering at co-rotation due to transient spi- 
ral density waves, which transfers stars from inner (or outer) 
disk regions into the solar vicinity by changing their angular 
momenta without alteration of their orbital circularity (hence 
eccentricities). These authors suggested that churning is the 
dominant process by which stars in the inner disk migrate out 
to the solar annulus, thus providing greater heterogeneity in 
the abundance and velocity distributions among solar neigh- 
borhood stars. 

The consequence of incomplete mixing from blurring and 
churning is that the metal-rich stars in the thin disk possess 
relatively lower rotational velocities (V^), while the metal- 
poor stars have higher V^. Thus, the expectation is that there 
should exist a trend of V<p with [Fe/H] among (at least) the 
thin-disk stars. Schonrich & Binney (2009a,b) indeed pre- 
dicted a significant downtrend of Vfy with [Fe/H] for the low 
[a/Fe] stars, due to incomplete mixing for younger stars. 
This prediction was confirmed by the later N-body models 
of Loebman et al. (2010), who employed slightly different 
treatments of radial mixing and star formation in their simu- 
lated disks from Schonrich & Binney (2009a,b), but found a 
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Table 1 

Summary of Observed Gradients for the Thin- and Thick-disk Subsamples 



Correlation 


AV /A[Fe/H] 


AV^/AR 


AVA|Z| 


Ae/A[Fe/H] 


Ae/AR 


Ae/A\Z\ 




(km s~' dex" 1 ) 


(km s kpc~' ) 


(km s kpc -1 ) 


(dex" 1 ) 


(kpc- 1 ) 


(kpc- 1 ) 


Thin 


-22.6±1.6 


-0.1 ±0.6 


-9.2 ± 1.2 


+0.000 ±0.005 


-0.003 ±0.002 


+0.036 ±0.003 


Thick 


+45.8 ±2.9 


-5.6 ± 1.1 


-9.4 ± 1.3 


-0.192 ±0.010 


+0.030 ±0.004 


+0.034 ± 0.004 



gradient of -19.7 km s" 1 dex" 1 for younger stars (identified 
with the thin-disk component with low [a/Fe]) in the solar 
neighborhood (7 < R < 9 kpc and 0.5 < \Z\ < 1 kpc). 

Our observed gradient of with [Fe/H] for the thin- 
disk stars in the range 0.5 < |Z| < 1.0 kpc (— 23. 5 ± 1.9 km 
s" 1 dex" 1 ) is not far from the estimate of -20 km s" 1 dex" 1 ob- 
tained by Loebman et al. (2010) for their simulated sample of 
young, low-[a/Fe] stars in their transition zone, which covers 
the same interval in height above the plane. Note that the scale 
of their [a/Fe] determinations and ours are slightly different, 
and they also employed the predicted oxygen abundance ratio 
as a proxy for [a/Fe], rather than the averages employed in 
our estimates. It appears that an overall velocity gradient of 
-20 to -30 km s" 1 dex"'with metallicity for the thin-disk sub- 
sample qualitatively agrees well with the expectations from 
the radial migration models. 

The extended tail of low- [a/Fe] metal-poor stars observed 
in Figures 1, 2, 4, and 5 can be also explained by the radial 
migration scenario. Roskar et al. (2008b) found that their 
simulated disk stars (when allowed to mix radially) exhibited 
a MDF more like the observations by Holmberg et al. (2007) 
than that obtained from an in situ sample without radial mi- 
gration. These authors concluded that radial migration was 
the likely cause of the broader MDF, which is also supported 
by our present data. 

It is noteworthy that our observed negative gradient of 
with [Fe/H] for the thin-disk stars (-20 to -30 km s" 1 dex" 1 ) 
stands in contradiction to expectations from traditional local 
evolution models in the solar neighborhood (without allow- 
ing for mixing or migration of stars), which predict a posi- 
tive slope of Vtj, with [Fe/H]. According to these models, the 
stars that were born early in the history of star formation in 
the thin disk are expected to be relatively metal-poor. These 
old metal-poor thin-disk stars should have experienced more 
perturbations, such as from variations in the Galactic poten- 
tial over time. As a result, such stars are expected to ex- 
hibit slower rotational velocities and larger velocity disper- 
sions than the younger, more metal-rich thin-disk stars. This 
inevitably leads to the expected production of a positive gra- 
dient of V<j, with [Fe/H], which we clearly do not find. 

When considering the radial migration models for the thick 
disk the case differs somewhat, in particular due to the much 
older ages of these stars. According to these models, it is ex- 
pected that most of the thick-disk stars that exist in the solar 
neighborhood today were born with high velocity dispersion 
in the inner portion of the Galaxy, in regions of higher local 
density, at a time when the metallicity of the ISM was rel- 
atively low and the a-abundance ratios were high. As they 
migrated outward over time, the lower gravitational restor- 
ing force of the local disk allowed these stars to explore or- 
bits reaching higher above or below the plane. Relatively few 
thick-disk stars are thought to have migrated inward from the 
outer disk region. These old stars had more time to experi- 
ence mixing of their orbits; in the case of complete mixing 
for these older stars, one might expect little or no trends be- 



tween rotational velocity and metallicity. 

Schonrich & Binney (2009a,b) did not make predictions 
of velocity trends with metallicity, on the grounds that in- 
sufficient knowledge of the earliest phases of disk formation 
exists to constrain expectations for such a potential gradient 
(i.e., unknown initial conditions). However, Loebman et al. 
(2010) reported from their simulation an insignificant gradient 
of+1.4kms _1 dex" 1 for these older stars (> 7 Gyr, which gen- 
erally matched the observed properties of the thick-disk com- 
ponent, e.g., high [a/Fe] ratios). Even though the migration 
strength in their simulation induced substantial mixing, the 
process was still incomplete. Thus, it would allow for the con- 
servation of significant velocity/metallicity trends. It should 
be mentioned, however, that their model was not specifically 
intended to match the properties of the Milky Way. 

In any case, the small or absent predicted correlations be- 
tween V<j, and [Fe/H] for the high- [a/Fe] stars from the mi- 
gration models are in contrast to our determination of a steep 
gradient of +40 to +50 km s" 1 dex" 1 for the observed high- 
ly a/Fe] stars we associate with the thick disk, as shown in Fig- 
ure [7] Thus, this trend of V<j, with [Fe/H] for the thick-disk 
subsample can provide a useful constraint to the radial migra- 
tion models mentioned above. 

In summary, the observed correlations between V$ and 
[Fe/H] for our low-[a/Fe] (thin-disk) stars can be naturally 
explained by the radial migration of stars from the outer disk 
(more metal-poor stars) and from the inner disk (more metal- 
rich stars) into the solar vicinity, as predicted by the migra- 
tion models. As explained by Schonrich & Binney (2009b), 
such a velocity gradient arises from the interplay between the 
churning and the blurring processes. The behavior of the 
high- [a/Fe] (thick-disk) stars is rather different (exhibiting 
a much steeper gradient) than expected from the simulated 
high-[a/Fe] stars of Loebman et al. (2010), although there 
remains the uncertainty of how well thick-disk stars are rep- 
resented in these models, and how well the models match the 
actual history of the Milky Way. It appears that stellar radial 
migration may have played an important role in the evolution 
of the thin disk, but, based on the information available from 
the current radial migration models and simulations, it is dif- 
ficult to ascertain the relative importance of radial migration 
for the formation and/or evolution of the thick disk. 

5.2. Rotational Velocity Gradients with Distance from the 
Galactic Center and Galactic Plane 

The gas-rich merger model of Brook et al. (2007) predicts a 
correlation between and R for stars in the disk system. Ac- 
cording to their simulations (especially their Figure 5) there 
should exist a detectable velocity gradient for the thin disk in 
the region of the solar neighborhood (R = 7-10 kpc). This dif- 
fers from our null gradient for the thin-disk subsample. Their 
simulations also indicate a negligible gradient for their thick- 
disk stars (which they refer to as "merger stars"), which is at 
least qualitatively in line with our small value of -5.6 ± 1.1 
km s" 1 kpc" 1 . However, it should be kept in mind that, as 
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Richard et al. (2010) demonstrated in their various gas-rich 
merger simulations, the initial orbital parameters of the merg- 
ers strongly affect the final kinematics and structures of the 
resulting disk populations. Brook et al. (2007) performed a 
simulation with a particular set of parameters to produce their 
disk systems, which may not necessarily match those of the 
Galaxy. For example, one rather large difference between this 
particular simulation and our results is that, while we find a 
difference in velocity lag of about 30 km s" 1 between our thin 
and thick-disk subsamples, the TV-body prediction calls for a 
difference of over 150 km s" 1 . Additional simulations of this 
process, better matched to the nature of the Milky Way, would 
clearly be useful to compare our results with. 

The dynamical heating of a pre-existing thin disk, as mod- 
eled by the simulations from Villalobos et al. (2010) also pre- 
dicts gradients of with respect to both R and |Z|. Looking 
at their Figure 14, the thickened-disk component exhibits a 
very weak trend of with R for low initial orbital inclina- 
tion of the merging satellite, while the correlation between 
the two quantities becomes stronger as the incidence angle 
is increased. Concerning the gradient of V<f, with \Z\, it is 
evident in their Figure 14 that the vertical gradient of is 
much shallower at higher orbital inclination. Thus, roughly 
speaking, our radial gradient of V<p for the thick-disk subsam- 
ple agrees better with that expected for low orbital inclina- 
tion of the merging satellite, but our vertical gradient is bet- 
ter matched by mergers with high orbital inclination. This 
may indicate that the compromise case of intermediate or- 
bital inclination (z = 30°) best describes our observed results, 
a possibility also considered by Villalobos et al. (2010). In 
any event, the comparisons of our thick-disk subsample with 
this particular model prediction imply that if the heating sce- 
nario played a major role in the formation of the thick disk, 
the initial orbital inclination of the merging satellite could not 
have been too small or too large. Of course, it is also possible 
that multiple satellite mergers may have been involved, which 
complicates these simple comparisons with a single merger. 

In summary, comparisons of our data with the gas-rich 
merger model from Brook et al. (2007) suggest that, while 
this model may not explain the lack of a rotational velocity 
gradient with Galactocentric distance for the thin disk, it does 
account for that observed for the thick disk. However, the 
much larger difference in the velocity lag than our finding be- 
tween the thin- and thick-disk stars remains to be resolved. 
The model of thin-disk heating by mergers of Villalobos et al. 
(2010) qualitatively agrees with the expected kinematic fea- 
tures of our thick-disk subsample, assuming that the merging 
satellite has an intermediate orbital inclination. 

In the previous section, which considered correlations be- 
tween Vtf, and [Fe/H], our results for the thin-disk population 
were shown to be in qualitative accord with predictions of the 
radial migration models, while those for our thick-disk popu- 
lation might not be. In order to be confident of the implica- 
tions of this result, one would like to compare with the predic- 
tions from more refined radial migration models that better re- 
produce the observed properties of the thick-disk population. 
On the other hand, the relationship between with R and 
|Z| for the high-[a/Fe] stars agrees better with the predictions 
of the gas-rich merger and thin-disk heating models that we 
have considered here. Taken as a whole, the presently avail- 
able comparisons of the various observed gradients suggest 
that the thick disk may have formed from either the mergers 
of gas-rich systems or the heating of a pre-existing thin disk 



by mergers, and has been little influenced by the secular pro- 
cess of stellar migration, while radial migration may well have 
strongly affected the evolution of the thin disk. 

5.3. Distribution of Stellar Orbital Eccentricities 

Sales et al. (2009) demonstrated that the orbital eccentrici- 
ties of a stellar population could also be used as a tool to probe 
the formation and evolution mechanisms of the disk system. 
In particular, taken at face value (and recognizing that their 
summary only pertains to a limited set of model parameters 
and histories), their Figure 3 suggests that radial migration 
models (e.g., Roskar et al. 2008a) generate symmetric dis- 
tributions of stellar eccentricities with rather narrow widths, 
while the gas-rich merger models (e.g., Brook et al. 2004, 
2005) produce distributions that are skewed toward higher 
eccentricity with larger widths. The accretion models (e.g., 
Abadi et al. 2003) distribute the eccentricities rather broadly 
over a wide range. For the disk heating scenario (e.g., Villalo- 
bos & Helmi 2008), there is a similarity of the eccentricity 
distribution with that of the merger model for e < 0.6, but 
there exists a secondary peak at high eccentricity (e ~ 0.8). 
Generally, they found that violent models such as disk heat- 
ing and accretion generated a distribution of stellar orbital ec- 
centricities spanning a large range, with secondary peaks at 
higher eccentricity, or at least with rather broad distributions 
of high eccentricity stars. By contrast, the smooth transition 
models, such as radial migration or in situ star formation from 
gas-rich mergers produced distributions dominated by lower 
eccentricity orbits covering relatively narrower ranges. 

Several studies have compared the above expectations from 
these models to observed distributions of orbital eccentricities 
for thick-disk stars in the solar neighborhood. Wilson et al. 
(201 1), for example, investigated the eccentricity distribution 
of a sample of thick-disk stars from RAVE. They concluded 
that their observed distribution, which peaked at low eccen- 
tricity and exhibited a lack of high eccentricity stars, disfa- 
vored the pure accretion model of Abadi et al. (2003), and 
was most consistent with the predictions of gas-rich merger 
models. Dierickx et al. (2010) carried out a similar test, us- 
ing a large sample of dwarfs from SDSS DR7, and suggested 
that their sample favored the gas-rich merger scenario as well. 
Casetti-Dinescu et al. (2011) performed an analysis using a 
sample of ^4400 red clump thick-disk stars from RAVE Data 
Release 2 (Zwitter et al. 2008) with available proper motions 
from SPM4. Their comparison of the derived orbital eccen- 
tricity distribution with model predictions supported the gas- 
rich merger scenario, or possibly the minor merger heating 
model (arguing that the expected secondary peak at high ec- 
centricity could be avoided, depending on the initial orbital 
configuration of the merging satellite(s)). Indeed, a recent 
simulation study by Di Matteo et al. (201 1) showed that, with 
the adoption of a particular set of initial conditions (a 1:10 
mass ratio and direct orbit of a presumed single interacting 
satellite), the disk heating model could also produce the dis- 
tribution of eccentricities observed by Wilson et al. (2011) 
and Dierickx et al. (2010) without creating a secondary peak 
at high eccentricity, confirming that the heating model may 
also be a viable mechanism for thick-disk formation. 

It is noteworthy that the various observational studies, 
based on different samples, with different distance estimates, 
and adopting different models for the Milky Way potential, 
all produce similar eccentricity distributions for the thick-disk 
population - a broad peak at low eccentricity and a lack of 
high eccentricity stars. Considering all of the studies men- 
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Figure 10. Left panel : Normalized distributions of eccentricities for the low-[a/Fe] stars (top left panel), high-[a/Fe] stars (top right panel), and for the entire 
G-dwarf sample (two bottom panels) over different |Z ranges. The vertical dashed line provides a reference at e = 0.2. The error bars are calculated by drawing 
from the sample with replacement 1000 times. A typical error in e is 0.044. Right panel : A reproduction of Figure 3 from Sales et al. (2009). The formation 
scenarios depicted are adopted from Abadi et al. (2003) for the accretion model, Villalobos & Helmi (2008) for the heating model, Roskar et al. (2008a) for the 
radial migration model, and Brook et al. (2004, 2005) for the gas-rich merger model. We adopt a bin size of 0.075 in the left-hand panels, close to that used 
in the Sales et al. figure. Taking the scale height (Zo) of the thick disk as 0.8 kpc, the range 1 < \Z/Zq\ < 3 in the Sales et al. figure corresponds to the range 
0.8 < |Z| < 2.4 kpc. 



tioned above, the favored mechanisms for thick-disk forma- 
tion are likely to be either (or both) the gas-rich mergers 
model or the thin-disk heating by minor mergers scenario, at 
least when considering only stellar orbital eccentricities as a 
probe. All of these studies rejected the pure accretion model 
of thick-disk formation (as advocated by Abadi et al. 2003). 

Unlike the previous observational studies mentioned above, 
which selected thick-disk stars mostly on the basis of spatial 
extent, we have selected a subsample of likely thick-disk stars 
based on their measured [a/Fe], as described in Section 3. We 
now compare the eccentricity distribution of our thick-disk 
subsample with expectations from each model cited in Sales 
et al. (2009). The left-hand set of panels of Figure [TOldisplavs 
the normalized distributions of eccentricity for the low- [a/Fe] 
(top left) and high-[a/Fe] (top right) populations. Each distri- 
bution in the top two panels is restricted to different slices on 
distance from the Galactic plane, as listed in the figure legend. 
The e distribution of the entire G-dwarf sample (without splits 
based on [a/Fe]), divided into regions that should emphasize 
the thin- and thick-disk regions, is shown in the bottom left 
and bottom right panels, respectively. 

The eccentricity distributions of the thin-disk subsample 
peak at much less than e = 0.2, with narrow widths, and ap- 
parently include very few high eccentricity stars (<? > 0.4) for 
the two |Z| regions shown in the top left panel. In contrast, 
the distributions for the thick-disk subsample shown in the 
top right panel peak at e ^0.2, and exhibit extended tails of 
higher eccentricities up to e ^0.8; there remains a relative 
lack of high eccentricity stars (e > 0.6). Although we find 
that the relative frequency of the high-e stars increases a bit 
at larger |Z| distance (red dashed line) for both subsamples, 
the distributions otherwise do not change significantly. This 
again confirms that the population separation based on [a/Fe] 
appears to work quite well. The eccentricity distributions for 
the full sample of G-dwarf stars exhibits some rather interest- 
ing features. Even at large |Z| distance (0.8-2.4 kpc, bottom 



right), where the thick-disk stars should dominate, the eccen- 
tricity does not appear similar to that of the thick-disk sub- 
sample separated by [a/Fe] (top right panel) in either range 
of |Z| distance; the peak and the width do not match. This 
underscores once more that, for the purpose of the selection 
of thick-disk stellar samples, purely spatial separations are in- 
sufficient. 

Comparing with the published model predictions in Sales 
et al. (2009), as shown in the right-hand panels of Figure 
[TOl the relative shortage of high eccentricity stars and the ab- 
sence of the secondary peak at high e ^0.8 in our observed 
distribution exclude the accretion origin and the disk heating 
model for the thick disk. Although the distribution expected 
from the radial migration models provides a viable descrip- 
tion of stars in the low eccentricity region, it fails to capture 
the observed high eccentricity tail of the thick-disk stars. The 
skewed distribution of observed eccentricities toward higher 
values is not well-represented by the radial migration predic- 
tions, which exhibit a more Gaussian-like shape (lower left 
panel of the right-hand panels of Figure [Tot . It should be 
noted, however, that an alternative radial migration model 
by Schtinrich & Binney (2009a,b) indicates the presence of 
a peak eccentricity between 0.1 and 0.2, with an extended 
tail towards high eccentricities, which is consistent with the 
shape of the observed e distribution of our thick-disk sample. 
Hence, we must be cautious in drawing firm conclusions on 
the formation mechanisms of the thick disk due to their appar- 
ent sensitivity to details of the models and simulations. Solely 
based on comparisons with the predictions in the published 
models from Sales et al. (2009), it seems that our eccentric- 
ity distribution most closely resembles that predicted from the 
gas-rich merger scenario. 

The eccentricity distribution of our thick-disk subsample 
differs little from the disk heating model of the Di Matteo 
et al. (2011) simulation. As identified by this simulation 
(and also mentioned in the discussions of the observations 
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Table 2 

Results of Qualitative Comparison Tests with Predictions of Published Models for Thick-disk Formation 



Model AV /A[Fe/H] A^/AS AV,<,/A|Z| e Distribution 

Accretion N/A N/A N/A Failed 

Disk Heating N/A Passed Passed Failed 

Radial Migration Indecisive 1 N/A N/A Indecisive 

Gas-rich Mergers N/A Passed N/A Passed 

Note. — N/A indicates that a model prediction is not available. The adopted thick-disk formation models are drawn from: Abadi et al. (2003) for the accretion model, Villalobos & 
Helmi (2008) for the heating model, Roskar et al. (2008a) for the radial migration model, and Brook et al. (2004, 2005) for the gas-rich mergers model. 

1 Based on the comparisons with predictions by Schonrich & Binney (2009a, 2009b) and Loebman et al. (2010). 



of Casetti-Dinescu et al. 2011 and Wilson et al. 2011), the 
secondary peak or high eccentricity region (e > 0.6) in the 
disk heating model is mostly occupied by the accreted stars 
(which retain the initial orbital characteristics of the merging 
satellite). In addition, depending on the initial conditions (es- 
pecially the orbital inclination of the interacting satellite) of 
the simulation, the high eccentricity secondary peak may (or 
may not) be seen in the predicted distribution of the eccentric- 
ities. In particular, the small satellite mass (1:10 mass ratio) 
in the Di Matteo et al. simulation would likely not contribute 
large numbers of stars to the solar neighborhood. 

It appears from our observed eccentricity distribution, and 
that of others, that the inclination of the merging small galaxy 
in the simulation of the disk heating model in Sales et al. 
(2009) may be less than 30°. This is qualitatively consistent 
with the findings from the correlations of V$ with R and |Z| in 
the previous section. We stress that the existence of the sec- 
ondary peak at high eccentricity and clear identification of the 
extended tail of the high eccentricity with observational data 
can provide strong constraints on the initial conditions on the 
merger or heating models. 

The observed eccentricity distribution of the full G-dwarf 
sample at larger distances from the Galactic plane, or for the 
thick-disk population with high [a/Fe] ratios, rule out a broad 
peak at intermediate eccentricity. This argues strongly against 
the importance of an accretion origin of the thick-disk com- 
ponent, unless the accretion model can explain the dominant 
population in the low eccentricity regime. At this stage, con- 
fident distinction between the other published models is in- 
feasible because uncertainties in the initial conditions, in the 
Af-body simulations themselves (e.g., artificial heating), and 
in the assumed model parameters (potential, secular heating, 
star formation histories, etc.) can produce differences that are 
roughly comparable to the predicted differences between var- 
ious scenarios. 

6. SUMMARY AND CONCLUSIONS 

We have assembled a sample of ^17,300 G-type dwarfs 
with available low-resolution (R ^2000) spectroscopy from 
SEGUE, a sub-survey conducted during SDSS-II. The sample 
we considered comprises stars with d < 3 kpc, log g > 4.2, 
spectra having S/N > 30 A -1 , V > +50 km s" 1 , [Fe/H] > 
—1.2, and 7 < R < 10 kpc. A separate test was carried out to 
eliminate a small number of stars that had larger probability 
of being associated with the halo than the disk system. 

Unlike the conventional assignment of stars into thin- and 
thick-disk components based on kinematics (or spatial distri- 
bution), we have made use of [a/Fe] as a reference to chem- 
ically divide our G-dwarf sample into likely thin-disk and 
thick-disk populations. 

Our chemically separated populations indicate that a nega- 
tive rotational velocity gradient with increasing [Fe/H] exists 



for the thin-disk population (-22.6 ±1.6 km s 1 dex *), while 
the thick-disk population exhibits a positive slope (+45.8 ± 
2.9 km s" 1 dex" 1 ) in the range 0.1 < |Z| < 3.0 kpc and R = 
7-10 kpc. Larger mean orbital radii are also noticed among 
the metal-poor thin-disk stars, as compared to the more metal- 
rich thin-disk stars, and smaller mean orbital radii are found 
for the thick-disk stars compared with the thin-disk stars. 

The distribution of rotational velocity appears independent 
of R for our thin-disk subsample, while there exists a very 
small correlation (-5. 6 ±1.1 kms" 1 kpc -1 ) between and/? 
for our thick-disk subsample. 

We have found that the observed lag of for the high- 
[a/Fe] stars relative to the low-[a/Fe] population is quite con- 
stant at a given |Z| distance (30 km s" 1 ), implying that our 
chemically separated populations are indeed distinct compo- 
nents in terms of their kinematics. This also allows us to infer 
that division by chemistry reveals the kinematic structure of 
each population better than division on the basis of spatial 
separation. 

The vertical gradient of with |Z| for our thick-disk sub- 
sample (-9.4 ± 1.3 km s" 1 kpc" 1 ) is smaller than that reported 
by Casetti-Dinescu et al. (2011; -25 km s" 1 kpc" 1 ), Ivezic 
et al. (2008; -29 km s" 1 kpc" 1 ), Girard et al. (2006; -30 
km s" 1 kpc" 1 ), and Chiba & Beers (2000; -30 km s" 1 kpc" 1 ). 
Without application of our proposed separation of the thin- 
and thick-disk subsamples, we find a vertical gradient of 
-19.4 ± 1.8 km s" 1 kpc" 1 for the stars with |Z| > 1.0 kpc, 
in better agreement with the previous observational studies. 
Hence, depending on how the thick-disk stars are selected 
(spatially or chemically), the kinematic trends change. As- 
suming the thick disk is a distinct component comprised of 
stars with high- [a/Fe] ratios, we again stress that the chemi- 
cal separation of the thick disk provides a more clear picture 
of the kinematics. 

It appears that there is no correlation between orbital ec- 
centricity and metallicity for the thin-disk subsample, while 
the trend of e for the thick-disk subsample rather steeply in- 
creases as the metallicity decreases. The e distribution for 
the low- [a/Fe] stars appears to be independent of R, whereas 
the high- [a/Fe] stars exhibit an increasing trend with distance 
from the Galactic center. The difference in average orbital 
eccentricity between the low- [a/Fe] and high- [a/Fe] subsam- 
ples appears constant at any given |Z| (about 0.1), which also 
indicates a clear distinction between these populations. Our 
approach of separating the thin-disk and thick-disk compo- 
nents by chemical tagging on their [a/Fe] abundance ratios 
yields well-defined, and distinct, kinematic trends for these 
populations as listed in Table Q] 

The rotational velocity gradient for the thin-disk subsam- 
ple with metallicity qualitatively agrees with the predictions 
of the radial migration models (Schonrich & Binney 2009b; 
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Loebman et al. 2010). Table [2] summarizes the results of 
qualitative comparisons of the various thick-disk formation 
scenarios with the observed properties of our G-dwarf sam- 
ple, based on predictions from these published (but still rather 
primitive) models. 

Based on these results, radial migration appears to have 
influenced the structural and chemical evolution of the thin 
disk, but may not have played a prominent role in the forma- 
tion and evolution of the thick disk. However, to be certain 
of this inference, comparisons with the predictions of more 
refined radial migration models that better reproduce the ob- 
served properties of the thick-disk population of the Milky 
Way galaxy are required. The preponderance of evidence, 
based on qualitative comparisons with existing thick-disk for- 
mation models, indicates that the thick disk of the Milky Way 
may have resulted from gas-rich mergers, or from heating of 
a pre-existing thin disk by minor mergers. We again empha- 
size that, although all of the models considered have had some 
success in reproducing aspects of the thick disk, no one the- 
ory has emerged that fully accounts for its detailed observed 
properties. We expect that newer generation models and sim- 
ulations will be strongly constrained by observations such as 
those presented here. 

Funding for SDSS-I and SDSS-II has been provided by 
the Alfred P. Sloan Foundation, the Participating Institutions, 
the National Science Foundation, the U.S. Department of En- 
ergy, the National Aeronautics and Space Administration, the 
Japanese Monbukagakusho, the Max Planck Society, and the 
Higher Education Funding Council for England. The SDSS 
Web Site is http ://www.sdss.org7| 

This work was supported in part by grants PHY 02-16783 
and PHY 08-22648: Physics Frontiers Center/Joint Institute 
for Nuclear Astrophysics (JINA), awarded by the U.S. Na- 
tional Science Foundation. J.A.J, acknowledges support from 
NSF grant AST-0607482. Z.I. acknowledges support from 
NSF grants AST 06-15991 and AST-07 07901, as weU as 
from grant AST 05-51161 to LSST for design and develop- 
ment activities. 

REFERENCES 

Abadi, M. G., Navarro, J. E, Steinmetz, M. & Eke, V. R. 2003, ApJ, 597, 21 
Abazajian, K., et al. 2009, ApJS, 182, 543 

Andrievsky, S. M., Luck, R. E., Martin, P., & Lepine, J. R. D. 2004, A&A, 
413, 159 

Aihara, H., et al. 2011, ApJS, 193, 29 
Allende Prieto, C, et al. 2008, AJ, 136, 2070 

An, D., Terndrap, D. M., Pinsonneault, M. H., Paulson, D. B., Hanson, R. 

B., & Stauffer, J. R. 2007, ApJ, 655, 233 
An, D., et al. 2009a, ApJ, 700, 523 
An, D., et al. 2009b, ApJ, 707, L64 
Barbanis, B., & Woltjer, L. 1967, ApJ, 150, 461 
Bensby, T., Feltzing, S., & Lundstrom, I. 2003, A&A, 410, 527 
Bensby, T., Feltzing, S., Lundstrom, I., & Ilyin, I. 2005, A&A, 433, 185 
Bond, N. A., et al. 2010, ApJ, 716, 1 

Brook, C. B., Kawata, D., Gibson, B. K., & Freeman, K. C. 2004, ApJ, 612, 
894 

Brook, C. B., Gibson, B. K., Martel, H. & Kawata, D. 2005, ApJ, 630, 298 
Brook, C, Richard, S., Kawata, D., Martel, H., & Gibson, B. K. 2007, ApJ, 
658, 60 

Carollo, D., et al. 2010, ApJ, 712, 692 

Casetti-Dinescu, D. I., Girard, T. M., Korchagin, V. I., & van Altena, W. F. 

2011, ApJ, 728,7 
Chiba, M., & Beers, T. C. 2000, AJ, 119, 2843 

Di Matteo, P., Lehnert, M. D., Qu, Y, & van Driel, W. 201 1, A&A, 525, L3 



Dierickx, M., Klement, R. J., Rix, H.-W., & Liu, C. 2010, ApJ, 725, L186 

Fuchs, B. 2001, MNRAS, 325, 1637 

Fuhrmann, K. 1998, A&A, 338, 161 

Fuhrmann, K. 2008, MNRAS, 384, 173 

Gilmore, G., & Reid, N. 1983, MNRAS, 202, 1025 

Gilmore, G., Wyse, R. F. G., & Bryn Jones, J. 1995, AJ, 109, 3 

Girard, T. M., Korchagin, V. I., Casetti-Dinescu, D. I., van Altena, W. F, 

Lopez, C. E., & Monet, D. G. 2006, AJ, 132, 1768 
Girard, T. M., et al. 2011, AJ, 142, 15 
Haywood, M. 2001, MNRAS, 325, 1365 
Haywood, M. 2008, MNRAS, 388, 1 175 

Holmberg, J., Nordstrom, B., & Andersen, J. 2007, A&A, 475, 5 19 
Ivezic, Z. et al. 2008, ApJ, 684, 287 
Jenkins, A. 1992, MNRAS, 257, 620 

Kazantzidis, S., Bullock, J. S., Zentner, A. R., Kravtsov, A. V., & 

Moustakas, L. A. 2008, ApJ, 688, 254 
Lemasle, B., Francois, P., Bono, G, Mottini, M., Primas, F., & Romaniello, 

M. 2007, A&A, 467, 283 
Lee, Y. S., et al. 2008a, AJ, 136, 2022 
Lee, Y. S., et al. 2008b, AJ, 136, 2050 
Lee, Y. S., et al. 2011, AJ, 141,90 

Loebman, S. R., Roskar, R., Debattista, V. P., Ivezic, Z., Quinn, T., & 

Wadsley, J. 2010, arXiv: 1009.5997 
Magrini, L., Sestito, P., Randich, S., & Galli, D. 2009, A&A, 494, 95 
Majewski, S. R. 1993, ARA&A, 31, 575 
Minchev, I., & Famaey, B. 2010, ApJ, 722, 112 
Munn, J., et al. 2004, AJ, 127, 3034 
Munn, J., et al. 2008, AJ, 136, 895 

Navarro, J. F, Abadi, M. G., Venn, K. A., & Freeman, K. C. 201 1, MNRAS, 
412, 1203 

Nissen, P. E., & Schuster, W. J. 2010, A&A, 511, L10 
Nordstrom, B., et al. 2004, A&A, 418, 989 

Parker, J. E., Humphreys, R. M., & Beers, T. C. 2004, AJ, 403, 74 
Quinn, P. J., Hernquist, L., & Fullagar, D. P. 1993, ApJ, 403, 74 
Reddy, B. E., Lambert, D. L., & Allende Prieto, C. 2006, MNRAS, 367, 
1329 

Reddy, B. E. 2010, in IAU Symp. 265, Chemical Abundances in the 

Universe: Connecting First Stars to Planets, ed. K. Cunha, M. Spite, & B. 

Barbuy (Cambridge: Cambridge Univ. Press), 289 
Richard, S., Brook, C. B., Martel, H, Kawata, D., Gibson, B. K., 

Sanchez-Blazquez, P. 2010, MNRAS, 402, 1489 
Robin, A. C, Reyle, C, Derriere, S., & Picaud, S. 2003, A&A, 409, 523 
Rocha-Pinto, H. J., Rangel, R. H. O., Porto de Mello, G. F, Braganca, G. A., 

& Maciel W. J. 2006, 453, L9 
Roskar, R., Debattista, V. P., Stinson, G. S., Quinn, T. R., Kaufmann, T., & 

Wadsley, J. 2008a, 675, L65 
Roskar, R., Debattista, V. P., Quinn, T. R., Stinson, G. S., & Wadsley, J. 

2008b, ApJ, 684, L79 
Sales, L. V., et al. 2009, MNRAS, 400, L61 

Schlesinger, K. J., Johnson, J. A., Lee, Y. S., Masseron, T., Yanny, B., 
Rockosi, C. M., Gaudi, B. S., & Beers, T. C. 2010, ApJ, 719, 996 

Schonrich, R., & Binney, J. 2009a, MNRAS, 396, 203 

Schonrich, R., & Binney, J. 2009b, MNRAS, 399, 1145 

Schonrich, R., Binney, J., & Dehnen, W. 2010, MNRAS, 403, 1829 

Sellwood, J. A., & Binney, J. J. 2002, MNRAS, 336, 785 

Sesar, B., Ivezic, Z., & Juric, M. 2008, ApJ, 689, 1244 

Sestito, P., Bragaglia, A., Randich, S., Pallavicini, R., Andrievsky, S. M., & 
Korotin, S. A. 2008, A&A, 488, 943 

Steinmetz, M., et al. 2006, AJ, 132, 1645 

Soubiran, C, Bienayme, O., & Siebert, A. 2003, A&A, 398, 141 
Soubiran, C, & Girard, P. 2005, A&A, 438, 139 

Spagna, A., Lattanzi, M. G., Re Fiorentin, P., & Smart, R. L. 2010, A&A, 
510, L4 

Spitzer, L., & Schwarzschild, M. 1953, ApJ, 118, 106 
Smolinski, J. P., et al. 2011, AJ, 141, 89 
Villalobos, A., & Helmi, A. 2008, MNRAS, 391, 1806 
Villalobos, A., Kazantzidis, S., & Helmi, A. 2010, ApJ, 718, 314 
Wilson, M., et al. 201 1, MNRAS, 413, 2235 

Wyse, R. F. G., Gilmore, G., Norris, J. E., Wilkinson, M. L, Kleyna, J. T., 

Koch, A., Evans, N. W., & Grebel, E. K. 2006, ApJ, 639, L13 
Yanny, B., et al. 2009, AJ, 137, 4377 
York, D. G., et al. 2000, AJ, 120, 1579 
Yoshii, Y. 1982, PASJ, 34, 365 
Zwitter, T., et al. 2008, AJ, 136, 421 



16 



Lee et al. 



APPENDIX 

RESOLVING OBSERVATIONAL CONFLICTS ON CORRELATIONS BETWEEN ROTATIONAL VELOCITY AND METALLICITY 

The recent study of Ivezic et al. (2008) concluded, on the basis of their analysis of a complete photometric sample from SDSS, 
that there existed little or no correlation between V<j, and [Fe/H] for stars in the disk system of the Milky Way, a finding confirmed 
by Bond et al. (2010). However, based on a spectroscopic sample of dwarfs from SDSS DR7, Spagna et al. (2010) reported 
a gradient of 40-50 km s" 1 dex -1 for stars with -1.0 < [Fe/H] < -0.5 and 1 < |Z| < 3 kpc. Loebman et al. (2010) claimed 
that the gradient found by Spagna et al. (2010) was caused by selection bias in the SDSS spectroscopic sample. As our current 
analysis also reveals a trend of with [Fe/H] for likely thick-disk stars, of similar size to that reported by Spagna et al., we have 
attempted to resolve these contradictory results between the various studies. 

As a first step, we employ the relationship devised by Ivezic et al. (2008) to obtain absolute magnitudes in the r band, and 
derive distances for our G-dwarf sample that should be on the same scale as theirs. Their adopted relationship is: 

M r (g-i, [Fe/H]) =-0.56+ 14.32*- 12.97 x 2 

+ 6.127x 3 -1.267x 4 + 0.0967x 5 (Al) 
-l.ll[Fe/H]-0.18[Fe/H] 2 , 

where x = (g — i). The above is the combined relationship of Equations A2, A3, and A7 in Ivezic et al. (2008). Then, by adopting 
the improved expression by Bond et al. (2010), we estimate the photometric metallicities for our sample. The adopted relationship 
is as follows: 

[Fe/H]pHOT =- 13.13+ 14.09x+28.04y-5.51xy 

-5.90x 2 -58.68y 2 + 9.14x 2 y (A2) 
-20.61xy 2 + 0.00x 3 + 58.20y 3 , 

where x = (u-g) and y = (g-r). All colors are reddening corrected, and note that the coefficient of the x 3 term is zero. We 
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Figure 11. Scatter plots (left panels) and histograms (right panels) of the differences in [Fe/H], distance, and V$, from top to bottom, between the photometric 
estimates (PHOT) and our estimates (SSPP) (see text for details of the determination of the "PHOT" parameters). 



Formation and Evolution of the Milky Way Disk System 



17 



S 
& 

H 
O 

K 

Oh 



240 

220 

S 200 h- 

| 180 

X 160 
> 

140 

240 
220 
200 
180 
~o 160 
140 

240 
220 
fl 200 



N = 3776 



-1.0 



1.0 < IZI SSPP < 1.5 kpc 



-*W 



-0.5 
[Fe/H] (SSPP) 



0.0 



Oh 
Oh 

00 



f N = 3938 




1.0 < IZI PHOT < 1.5 kpc ; 


* * 


■ 1 T 




-1.0 




-0.5 0.0 
[Fe/H] (PHOT) 


f N = 3781 




1.0 < IZI SSPP < 1.5 kpc ~ 


r * * 






-1.0 




-0.5 0.0 
[Fe/H] (PHOT) 



180 
160 
140 



Figure 12. Derived rotational velocity gradients as a function of metallicity. The top panel shows the trend with the SSPP-derived quantities for all stars at 
|Z| = 1.0-1.5 kpc from our G-dwarf sample. Note that the SSPP-determined \Z\ distance is used to select these stars. The middle panel are the photometric 
determin ation s in the same interval of height above the plane. Notice that there are more objects selected in this |Z| distance range, which was calculated using 
Equation lAll The bottom panel displays the trend of our values of vs. the photometric metallicity, and a very similar pattern between the middle and the 
bottom panel is observed. The error bar on each point is calculated by resampling 100 stars with replacement 1000 times. 



refer to the distance determined with Equation lAll as the "photometric distance", the metallicity estimated by Equation lA2l as the 
"photometric metallicity", and the rotation velocity calculated using the photometric distance in combination with the measured 
radial velocities and proper motions as the "photometric rotational velocity". The label "PHOT" in Figures [TTI and [121 indicate 
these estimates, while the label "SSPP" denotes the values we have used for the G-dwarf sample. 

FigureQT]shows scatter plots (left panels) and histograms (right panels) of the differences in [Fe/H], distance, and between 
the photometric (PHOT) estimates and our estimates (SSPP) for stars in |Z| = 1.0-1.5 kpc from our G-dwarf sample considered 
in this paper. It can be noticed from inspection of the top panels that the photometric metallicity is consistently higher at low 
[Fe/H] and lower at high [Fe/H], compared to the SSPP estimates, with an overall shift of about 0.1 dex (top right panel). 

The middle panels suggest that our distance determination is on average ~0. 1 kpc lower than the photometric estimate, with a 
trend that this deviation becomes larger as the distance increases. The bottom panels show that our rotational velocities generally 
agree with the photometric rotational velocities, but our values trend higher at low V$. 

With these differences (especially for [Fe/H]) kept in mind, we have examined the trend of V$ with [Fe/H] based on the 
photometric estimates and based on the spectroscopic estimates, for stars in the range |Z| = 1.0-1.5 kpc. Figure [121 shows the 
results. The top panel makes use of the spectroscopically-derived quantities, while the middle panel comes from the photometric 
quantities. Due to the small difference in the distance estimates (larger distances for the photometric distance), we see in the 
middle panel that there are more objects selected in this |Z| distance range, which does not affect our conclusions. It is obvious 
that we obtain a flattening of the V<f, relationship below [Fe/H] < -0.5 for the photometrically-determined values. Even if we 
consider the photometric metallicity and our values of Vj,, we obtain a very similar pattern (bottom panel). This strongly suggests 
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that the effect of the input metallicity in Equation [AT] on distance for calculation of is minimal. However, it is clear that the 
difference between use of the photometric metallicity and the spectroscopic metallicity makes a large difference in the derived 
trend of V^. This can be accounted for by the scattering of higher metallicity stars (which have high V$) into the photometrically 
determined low-metallicity region, resulting in a flatter gradient of V<j, with [Fe/H] for stars with [Fe/H] < -0.5, as seen by 
comparing the middle panel of Figure [12] with the top panel. 

There are at least two reasons that the photometric metallicity relation in Equati on IA2I may assign the metallicity of a star 
to a value that strongly deviates from that estimated from the SSPP. First, Equation I A2I was not calibrated with the metallicity 
estimates used in this study, but with those available in DR7. There has been significant improvement in the SSPP for estimation 
of [Fe/H] since the DR7 release, especially for stars with near-solar and super-solar values (see Smolinski et al. 2011). For a 
proper comparison, the photometric metallicity relation needs to be re-calibrated with the metallicities available from the DR8 
release. The other reason is that small random errors in the photometric measurements (and zero points) can strongly influence 
the photometric metallicity estimate (and its errors) (Z. Ivezic et al. in preparation). As a result, stars that are in reality of high 
metallicity can be artificially moved into the low-metallicity region. 



